Quantum phase transitions in superconducting arrays under external magnetic fields 
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We study the zero-temperature phase transitions of two-dimensional superconducting arrays 
with both the self- and the junction capacitances in the presence of external magnetic fields. We 
consider two kinds of excitations from the Mott insulating phase: charge-dipole excitations and 
single-charge excitations, and apply the second-order perturbation theory to find their energies. 
The resulting phase boundaries are found to depend strongly on the magnetic frustration, which 
measures the commensurate-incommensurate effects. Comparison of the obtained values with those 
in recent experiment suggests the possibility that the superconductor-insulator transition observed in 
experiment may not be of the Berezinskii-Kosterlitz-Thouless type. The system is also transformed 
to a classical three-dimensional XY model with the magnetic field in the time-direction; this allows 
the analogy to bulk superconductors, revealing the nature of the phase transitions. 

PACS numbers: 74.50.+r, 67.40.Db 



I. INTRODUCTION 

Two-dimensional (2D) superconducting arrays with 
charging effects have drawn much interest because of 
their interesting phase transitions jj| and relations to 
other systems such as the Bose-Hubbard model and the 
quantum XX Z spin model ||. In those arrays where 
both the self-capacitance Cq and the junction capaci- 
tance C\ are present, charging energies due to both ca- 
pacitances need to be considered simultaneously since 
the nature of the phase transition in general depends 
on them. For example, at zero temperature the system 
with only self-capacitance can be mapped into a classical 
three-dimensional (3D) XY model with the ratio Eq/Ej 
of the charging energy E = e 2 /2Co to the Josephson 
coupling strength Ej taking the role of the tempera- 
ture. In the opposite limit, on the other hand, it is well 
known that there exists an interesting duality between 
charges and vortices and the system with only junc- 
tion capacitance at zero temperature undergoes a charge- 
unbinding Berezinskii-Kosterlitz-Thoulcss (BKT) transi- 
tion [0 from the insulating phase to the superconducting 
one as the ratio Ej/E% with E\ = e 2 /2Ci is increased. 
The critical value (Ej/Ei) c , beyond which the array is 
superconducting, has been found 0.6 in experiment 
0.23 in the duality argument [Q and in the variational 
method M, 0.51 in perturbation expansion S, and 0.36 
in quantum Monte Carlo simulations j^]. As Co is in- 
creased from zero in this system, the interactions between 
charges are screened, with the screening length given by 
A = y/Ci/Co 0, and the nature of the phase transition is 
expected to alter. Recently, the zero-temperature phase 
diagrams have been studied by means of the perturbation 
expansion ||, which suggests that the phase transition 
from the Mott insulating phase to the superconducting 
phase is governed by the single-charge (SC) excitations 
or by the charge-dipole (CD) excitations, depending on 



the ratio Ci/Cq as well as on the charge frustration. 

In this paper, we study the quantum phase transi- 
tions in 2D superconducting arrays, focusing on the ef- 
fects of external magnetic fields as well as the compe- 
tition between self- and junction capacitances. The ef- 
fects of external magnetic fields on phase transitions have 
been studied both in the classical arrays without charg- 
ing energy |lj and in the quantum arrays |9|~p"2||. How- 
ever, most existing analytical works on the latter have 
employed mean-field-like approximations, which are not 
reliable in two dimensions. We thus adopt the pertur- 
bative expansion instead, and consider the SC and the 
CD excitations to the second order in Ej. The obtained 
phase boundaries between the Mott insulating phase and 
the superconducting one exhibit strong commensurabil- 
ity effects due to the magnetic frustration. The results 
are compared with experiments, suggesting the possibil- 
ity that the experimentally observed transition may not 
be of the charge-unbinding BKT type. It is also shown 
that the dual transformation maps the system with both 
magnetic frustration and general capacitance onto a clas- 
sical 3D XY model under the magnetic field in the time- 
direction. This transformation allows us to discuss the 
nature of the phase transitions, by analogy with bulk su- 
perconductors under magnetic fields. 

There are four sections in this paper: Section II is 
devoted to the perturbative approach, from which the 
zero-temperature diagrams are obtained. We compare 
the phase diagrams with those obtained from the mean- 
field approach and those observed in experiments. In 
Sec. Ill the system is transformed into a classical 3D 
XY model under the magnetic field in the (imaginary) 
time direction. The nature of the phase transitions is 
discussed by analogy with bulk superconductors under 
magnetic fields. Finally, Sec. IV summarizes the results 
and presents some discussions. 



II. PERTURB ATIVE APPROACH 

We begin with the Hamiltonian describing the super- 
conducting array with magnetic frustration: 



/ 

H Q + V, 



Ej cos ( 



A lj ) (1) 



where the charge qj at site j represented in units of 2e 
(e > 0) and the phase ^ of the superconducting or- 
der parameter at site k satisfy the commutation relation 
[<^fc)<Zj] = iSjk- Whereas the summation in the second 
term runs over all nearest neighboring pairs, the charges 
interact via the inverse of the dimcnsionless capacitance 
matrix defined by 



Cij = (1 + 4Ci/C )5 it j - (Ci/C )(5i ij+ ^ + <5ij-i + ^i,j+y + $id-y) 



(2) 



with the self-capacitance Co and the junction capacitance 
C\. The important external parameters E = e 2 /Co 
and Ej measure the self-charging energy and the Joseph- 
son coupling energy, respectively. The magnetic bond 
angle A^j between the two sites i and j is given by 
the line integral of the magnetic vector potential A: 
Aij = (27r/$o) J° A • dl with the magnetic flux quantum 
$ = hc/2e. 

When Ej = 0, the system described by the unper- 
turbed Hamiltonian Hq has the Mott insulating phase as 
the ground state with the charge configuration % = 
at any site. In the opposite limit of Eo = 0, on the 
other hand, the system is described by the 2D classical 
XY model, displaying superconductivity at zero temper- 
ature. It is thus expected that there exists a critical value 
of Ej j Eq beyond which the superconducting phase be- 
comes energetically favorable. The critical value may be 
determined by comparing the energy of the Mott insulat- 
ing phase with that of the superconducting phase. The 
energy Em of the Mott insulating phase is easily com- 
puted up to the second order in Ej/Eq: 



Em — —- 



EjN 



8E (C m - C A Q ) 



(3) 



for the L x L square array (N = L 2 ) (see Ref. Q for the 
detailed calculation). It is of interest to note here that 



Em does not depend on the external magnetic field. 

Since the lowest excitation that can lead to the change 
of the ground state from the Mott insulating into the 
superconducting phase is presumably point-like, and we 
consider two kinds of excitations: the SC and the CD 
excitations. We first consider the SC-type excited state 
and write its energy up to the second order in Ej/Eq: 



Eh 



sc 



E. 



(o) 
sc 



E. 



(i) 
sc 



E, 



(2) 

sc- 



(4) 



In the SC-type excited state a single positive charge is 
located only at one site, and the zeroth-order energy is 
given by 



Ef^ = ^Y^q^q^ = AE C 



/ 



(5) 



Since the single charge can be located at any site with- 
out any energy difference, we need to make use of the 
degenerate perturbation theory. In the first order it is 
thus necessary to diagonalize the matrix 



<i|V|j) = -{i\Ej 



E 

(k,i) 



$3 

co8(^> k -(t>i-A k i)\j) = — ^-Pij, (6) 



where \i) denotes the charge eigenstate with the single 
charge at site i and the matrix element Pij is defined by 



= J exp(-iAij) 
*ij — 



for nearest neighboring sites i and j, 
otherwise. 



In the presence of the external magnetic field, the gauge- 
invariant magnetic frustration is defined by / = <!>/ 4>o 
with $ being the magnetic flux per plaquette. If / = p/q 
with p and q being relative primes, it is obvious that 
is invariant under the magnetic translation of q lattice 
sites. Noting this translational symmetry, we represent 
the position of site i by the vector R+a, where R is the 



position vector of the q x q superlattice unit cell con- 
taining site i and a denotes the relative position of site i 
inside the superlattice (see Fig. |l|), and write 



P^ = P(R, a; R', a') = P(R-R, a; 0, a'). 



Through the Fourier transformation 



(7) 



2 



u( p )(a) = -^Y e 4p - R i;(R+a), 
^ R 



p(p)(a;a')^^e l P R P(R,a;0,a') 

R 



(8) 



the matrix is block-diagonalized, resulting in the eigen- 
value equations 



E 



P( p >(a;a> (p) (a') = A< p V p )(a), 



(9) 



where M = N/q 2 is the total number of superlattices, 
and v (R+a) = Vi is the wavefunction for the eigenstate 
of Rij. The numerical diagonalization of M Hermitian 
matrices p( p )(a;a'), each of the size q 2 x q 2 , yields the 
largest eigenvalue P max , which in turn gives the first- 
order energy 



-^SC — PjPmax/2. 



(10) 



It is well known that the largest eigenvalue P max always 
shows up for p = 0, and further, the eigenstates corre- 
sponding to the eigenvalue P max are g-fold degenerate 
for / = p/q p3[ . Since the degeneracy is not completely 
removed in the first-order calculation, it is necessary to 
diagonalizc the second-order matrix Q with the clement 
given by 



x ^ (rf|y|q)(q|FK) 
^ dd ' - 2^ p (0) _ p,(0) 

q^D -^sc c q 



(11) 



where the summation runs over all the charge eigenstates 
|q) outside the space D spanned by the SC states, and \d) 
represents the c?th eigenstate corresponding to the eigen- 



value Prr 



M>=E 



Vd,i\V 



E 

R,a 



v d (R+a) | R+a) 



(12) 



Here the wavefunction t>d(R+a) of the dth degenerate 
state is related to the component u,j(a) of the normal- 
ized eigenvector of P via 



v d (R+n) = 



<M 



(13) 



where the superscript in u p_0 (a) has been omitted for 
simplicity [see Eq. (fa)]. Inserting Eqs. ([IJ) and (|T^ ) into 
Eq. (p"l|), we obtain 



Qdd> = ^^<(a)v rf /(a')QR-i 



(14) 



R 



with 



VR+a,a' - 2^ (0) „(0) 



It is obvious that Qij does not vanish only when \i) 
and \j) are related by two successive charge hoppings 
as shown in Fig. |[ When i = j (denoted by an empty 
circle in Fig. 0), it is easy to find that 



1 



32E ° m ( c n? ~ - (Coo 1 - Cr ) 

where the summation y~)r m „) runs over m and its four 
nearest neighbors n with the restriction of nonzero de- 
nominator. When j = i ± x ± y (denoted by 'x' symbols 
in Fig. |^), there exist two intervening states |q), yielding 



K, (15) 



li-,3— z±x±y 



E 2 e i(A j ,i±*+A i ±5t,i) _j_ e i(A3,i±^A i±94 ) 

32P (C^-C^-iC^-C^) • 



(16) 



On the other hand, when j = i ± 2x(2y) as represented 
by filled circles in Fig. |2| the matrix element is computed 
to be 



Q 



ij=i±2x(2y) 



o*(^3',i±4(y)+^i±il(y),i) 



32^o (C^ -C-')- (Coo 1 -O 



(17) 



Equations (|lj) - (0) together with the eigenvector u 
obtained in the first-order calculation yield the explicit 
form of the matrix element Qdd' in Eq. (|l4]). The matrix 
Q is then diagonalized to give the minimum eigenvalue 
Qmini which in turn leads to the energy of the SC state: 



sc 



E, 



(o) 

SC 



^sc + ^sc = ^EqCqq — PjP max /2 + Q n 



? {2) 



(18) 



Comparing it with the energy of the Mott phase in 
Eq. (||), we find the phase boundary between the Mott 
insulating phase and the superconducting phase. Fig- 
ure H shows the obtained phase boundaries separating 
the superconducting phase (above each curve) and the 



insulating one (below) on the plane of Ej/Eq and /, for 
various values of C\/Cq. In obtaining the critical val- 
ues (Ej/E ) c numerically, we have considered systems of 
sufficiently large sizes, so that the critical values display 
convergence in at least three significant digits, for given 



3 



parameters / and C±/Cq. Thus the system size has been 
increased up to N — 384, depending on the values of / 
and Ci/Cq, and the convergence has been confirmed. 

It is obvious that the superconducting region expands 
as the junction capacitance C± is increased, confirming 
the results in Refs. |7]||. In particular, the obtained 
phase diagrams are entirely similar to those obtained for 
Co ^> C\ in the mean- field approximation [f^| , demon- 
strating significant commensurate-incommensurate ef- 
fects due to the magnetic frustration. Quantitatively, 
however, there exists discrepancy between the results of 
the perturbation expansion and those from the mean- 
field approaches: The estimated critical values in the 
former are rather larger. Furthermore, the perturba- 
tion expansion yields the ratio of the critical value for 
/ = 1/2 to that for / = approximately given by 1.9 for 
Ci/Co ^ 0.1 (see Fig. 3); this is larger than the value \/2 
predicted in the self-charging limit within the mean-field 
approximation The increase of C\/Cq is found 

to reduce the ratio monotonically. It is of interest to no- 
tice here that the first-order calculation reproduces the 
mean- field value \/2 regardless of C\/Cq. 

We now consider the CD-type excited state, where 
there exists a pair of positive and negative charges sepa- 
rated by the lattice spacing. The energy of the CD state 



is written as 

^cd = £ C ° D + £od + ^CD ( 19 ) 

up to the second order in Ej/Eq. The zcroth- and 
the first-order energies are easily calculated: i%D = 
8i?o(C 00 1 — C^l) and E^ = 0. To calculate the second- 
order term, we apply the second-order degenerate pertur- 
bation theory, and diagonalize the matrix M , the element 
of which is given by 

(i, 3 \M\k,l) = }^ (o) _ (o) • ( 2 °) 

q<£D ^CD 

Here \k,l) is the CD states with the positive charge at 
site k and the negative charge at site where / is one 
of the four nearest neighbors of k, and the sum is per- 
formed over the intervening states |q) outside the space 
D spanned by all CD states. The above matrix clement 
does not vanish only when \k, I) can be connected to 
by two successive charge hoppings. Figure [| shows all 
possible states of \k, I) when \i,j) is given as in Fig. 0(a). 
While the matrix element corresponding to Fig. |](a) is 
given by 



F 2 i 1 
(i i\M\i i) = J \ - 

* 2E »fa^-c^)-{cJ-c£)-{c^-c^y 

there is only one intervening state in the case of Fig. |J(b), leading to 



(21) 



(i,j\M\j,i) 



32£ C V - Cll 



For \k,l) given by the state in Figs. |J(c) and (d), we find 



(22) 



(i,j\M\i,l=i a ) 



32E Q 



s~i — 1 s « — 1 s~Y—l — 1 o ^ ' — 1 > ) s i — 1 s~i — 1 

o 00 ~ °x,0 ^x+y.O - U i,0 aO x,0 - ZLy x+y,0 - ^00 



^±,0 ~ ^x+y,0 ~~ 2( ^00 J 



(23) 



where i a denotes the crth nearest neighboring site of i (a = 1,2,3,4). Similarly, we get 



{i,j\M\k=j a ,j) 



32E 



■AAj,i+A kJ ) pi(Ai a ,i + A k , in ) 



oi(A h , ia +Ai a ,i) 



°00 °x,0 U x+y,0 U x,0 ,3L/ x,0 ~ ZO x+y,0 - O 00 

e i(A htj +A jii ) 



for Figs. ||(e) and (f), 



(24) 



(i ' i|M|M)= 32£b 



i(A s ,i+Au,i) 2i{A kii +A jt i) 



,i(A kt i+AjA) 



L°00 ~~ °x,0 ^x+y.O _ U x,0 ZO x+y,0 ~~ dO x,0 °00 



(25) 
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for Figs. 0(g) and (h), and finally 
(i,j\M\k,l) = 



El 



32En 



J{Aj,i+A ki i) 



°00 



J{Aj,i+A k ^i) 



(26) 



for the cases corresponding to Fig. [|(i). Equations (|2l|)- 
( p6| ) give the ANx4N matrix M, which, again via the 
Fourier transformation, reduces to Aq 2 x4q 2 Hermitian 
matrices for / = p/q. By diagonalizing numerically the 
resulting matrices, we obtain the second-order energy 

(2) 

-Eqq, the comparison of which with the energy of the 
Mott insulating phase given by Eq. (J||) yields the phase 
boundaries. Figure [5] displays the obtained boundaries 
in the plane of Ej/E[ and /. As in obtaining Fig. 3, we 
have considered systems of such large sizes that at least 
three significant digits of the critical values (Ej/Ei) c are 
obtained for given values of / and Ci /Co- 
in experiment on the arrays with Ci/Cq ~ 100, the 
critical values (Ej/Ei) c w 0.6 for / = and 0.9 for 
/ = 1/2 have been observed |(|; this is to be compared 
with the corresponding values obtained in the pertur- 
bation scheme here, with the CD excitations taken into 
consideration: {E 3 /E x ) c 0.503 for / = and 0.508 for 
/ = 1/2. Remarkably, the consideration of CD excita- 
tions yields the critical value for / = 1/2 not far larger 
than that for / = 0, in disagreement with the experi- 
mental result. On the other hand, the phase boundary 
computed from the consideration of the SC excitations 
for Cx/Cq = 100, shown in Fig. ^, is in general consistent 
with that measured in experiment Contrary to the 
usual anticipation that CD excitations play a crucial role 
in destroying the Mott insulating phase for C\ ^ Co, this 
apparently suggests that the superconductor-insulator 
transition observed in experiment is driven by SC ex- 
citations rather than by CD ones, raising the interesting 
possibility that the transition may not be of the BKT 



type. 



III. DUAL TRANSFORMATION APPROACH 

In this section, we examine the nature of the phase 
transitions discussed above in terms of the topological 
excitations. This is achieved by means of the mapping 
of the quantum phase model given by Eq. ([l]) into an 
effective 3D classical model; this approach was already 
adopted by other authors in the absence of the exter- 
nal magnetic field [Q. We begin with the Euclidean ac- 
tion, corresponding to the Hamiltonian in Eq. (|l|) , in the 
imaginary-time path-integral representation: 



(27) 



where K = ^Ej/8E , A^.j = Ajj +fl , V M (fi = x,y) 
and V r are difference operators in the space and time 
directions, respectively, and the (imaginary) time-slice 
spacing St has been chosen to be ^/8EjEq/?i In the 
nearest- neighbor charging limit (Co = 0), the coupling 
constant and the time-slice spacing in Eq. ( p7| ) are given 
by K = \J Ej /8Ei and St — \/8EjE[/h, respectively. 
Standard procedures (l^Jl^l then lead to the dual form 
of Eq. ( p7j ) , which is simply the effective Hamiltonian for 
the 3D classical system of vortex loops: 



H v = -ttI<Y^ t) - (S^r/l ^(r, - r,, t - t') [^(r 3 ,r') - ^ :T /| . 

i,j,r,r' fi=X,y,T 

Here the vortex interaction [/^(r, r) = 2ir [L^(0, 0) — C/ M (r, r)] is determined by the Fourier transforms 

C(q) 



(28) 



U x (q,w) = U y {q,u) 



U T (q,oj) 



AW + Ay+C(q)A(«) 

1 



A(q x ) + A(q y ) + C(q)AH 



(29) 
(30) 



with A (q) = 2(1— cosq). The Fourier transform of the ca- 
pacitance matrix is given by C(q) = 1 + (C\ / Co)[A(q x ) + 
A(q y )} for C ^ and C(q) = A(q x ) + A{q y ) for C = 0. 



Note also that the vortex lines can terminate nowhere 
but form closed loops or go to infinity, as implied by the 
condition V • v(r, r) = 0. 
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The behavior of the interaction C/ M (r, r) in Eq. 
depends crucially on whether Cq vanishes, since the 
Coulomb interaction between charges (Cooper pairs) 
is infinitely long-ranged for Co = 0. If Co ^ 0, 
one can show, in the same manner as in Ref. p6[ , 
that at large scales (V r 2 + t 2 3> 1) the interaction 
Up, is isotropic and displays the asymptotic behavior 
Ufj,(r, t) ~ —1/ yr* + t 2 apart from an additive constant, 
regardless of the ratio C\/Cq. Accordingly, the system 
is described by the 3D isotropic XY model under an ex- 
ternal magnetic field in the r-direction, the topological 
representation of which is given by Eq. J2g| ) . The 3D XY 
model has been widely used as a model for the bulk su- 
perconductor at temperatures low enough to ne glec t the 
amplitude fluctuations of the order parameter JL7J. By 
analogy with the vortex lattice melting transition at the 
temperature T m (H) in the mixed state of a type-II su- 
perconductor, a first-order phase transition is expected 
at K c (f) in our system under the magnetic field, as K is 
increased from zero p7| , p^| . At zero field, in particular, 
the phase transition should be continuous, belonging to 
the 3D XY universality class Q] . 

For Co = 0, on the other hand, the interaction is highly 
anisotropic: U x {y,t) — U y (r,T) ~ exp(— V r 2 + r 2 ) 
while U T (r, r) ~ e~' T ' logr. Thus the equivalent classical 
system described by Eq. (|2|) forms a layered structure of 
planar spins with strongly anisotropic coupling constants. 
The 3D anisotropic XY model has also been studied ex- 
tensively as a special case of the Lawrence-Doniach model 
for high-temperature superconductors |lj],|l9| : The effec- 
tively low dimensionality enhances the phase fluctuations 
and lowers the transition point K c (f ) Jl7[ ] . Furthermore, 
at zero field, the strong anisotropy drives the transition 
to be of the BKT type 

These arguments have been summarized in Fig. |t]. 
Note that the important effects of frustration arising from 
the commensurability between the flux lattice and the 
underlying lattice are disregarded here. Such continuum 
approximation is believed to be qualitatively valid in low 
field regions (represented by the solid lines in Fig. 0). 
As the field is increased, frustration effects are expected 
to come into play and to yield sensitive dependence of 
the transition on the field, reproducing the perturbative 
results shown in Figs. || and [| 



the commensurate-incommensurate effects. In particu- 
lar, the superconductor-insulator transition observed in 
experiment has been suggested to be driven by single- 
charge excitations rather than by charge-dipole ones, 
and thus the possibility that the transition may not 
be of the Berezinskii-Kosterlitz-Thouless type has been 
pointed out. In this regard, it is noteworthy that the low- 
est excitation has already been shown to be comprised by 
the single-charge type rather than the charge-dipole type 
even for large values of Ci/C , so long as there exists 
finite charge frustration Q. In experiment, it is practi- 
cally impossible to set the charge frustration exactly zero, 
and accordingly, the lowest excitation should presumably 
be of the single-charge type even in the nearest-neighbor 
charging limit. Indeed it has recently been pointed out 
that the absence of the Berezinskii-Kosterlitz-Thouless 
charge-unbinding transition in experiments |j] may be 
attributed to the presence of the finite charge frustration 
which is randomly distributed over the arrays |2(i(| . 

We have also transformed the system to a 3D classical 
XY model under magnetic fields in the time direction. 
The nature of the phase transitions at low magnetic fields 
has been discussed, based on the analogy with the bulk 
superconductor under magnetic fields. Unfortunately, 
unlike in the 2D XY model, there have been few stud- 
ies of the frustration effects in the 3D XY model, which 
disallows quantitative comparison at this stage. Never- 
theless, the analogy with the continuum superconductor 
provides a complement to the perturbative estimate of 
the phase boundaries, already giving a clue to the nature 
of the phase transitions. 
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IV. CONCLUSION 

We have studied the zero-temperature phase transi- 
tions of two-dimensional superconducting arrays with 
both self- and junction capacitances in the presence 
of external magnetic fields. Through the use of the 
second-order perturbation theory, we have considered 
both single-charge excitations and charge-dipole excita- 
tions, from which the phase diagrams are obtained. It 
has been found that the phase boundaries are quite sen- 
sitive to the variation of the magnetic frustration due to 



Present address: Department of Theoretical Physics, 
Umea University, S-901 87 Umea, Sweden 

[1] For a list of references, see, e.g., Proceedings of the 2nd 
CTP workshop on Statistical Physics: KT Transition and 
Superconducting Arrays, edited by D. Kim, J.S. Chung, 
and M.Y. Choi (Min-Eum Sa, Seoul, 1993); Physica B 
222 253-406 (1996). 

[2] A. van Otterlo, K.-H. Wagenblast, R. Baltin, C. Bruder, 
R. Fazio, and G. Schon, Phys. Rev. B 52, 16 176 (1995). 

[3] L.V. Keldysh, Pis'ma Zh. Eksp. Teor. Fiz. 29, 716 (1979) 
[JETP Lett. 29, 658 (1979)]; J.E. Mooij, B.J. van Wees, 



G 



L.J. Geerligs, M. Peters, R. Fazio, and G. Schon, Phys. 

Rev. Lett. 65, 645 (1990). 
[4] R. Fazio and G. Schon, Phys. Rev. B 43, 5307 (1991), 

and references therein. 
[5] V.L. Berezinskii, Zh. Eksp. Teor. Fiz. 59, 907 (1970) [Sov. 

Phys. JETP 32, 493 (1971)]; J.M. Kosterlitz and D.J. 

Thouless, J. Phys. C 6, 1181 (1973); J.M. Kosterlitz, ibid. 

7, 1046 (1974); J.V. Jose, L.P. Kadanoff, S. Kirkpatrick, 

and D.R. Nelson, Phys. Rev. B 16, 1217 (1977). 
[6] H.S.J, van der Zant, L.J. Geerligs, and J.E. Mooij, Eu- 

rophys. Lett. 19, 541 (1992); H.S.J, van der Zant, W.J. 

Elion, L.J. Geerligs, and J.E. Mooij, Phys. Rev. B 54, 10 

081 (1996). 

[7] B.J. Kim and M.Y. Choi, Phys. Rev. B 52, 3624 (1995). 
[8] B.J. Kim, J. Kim, S.Y. Park, and M.Y. Choi, Phys. Rev. 

B 56, 395 (1997). 
[9] C. Rojas and J.V. Jose, Phys. Rev. B 54, 12 361 (1996). 
[10] A. Kampf and G. Schon, Phys. Rev. B 37, 5954 (1988); 
A. van Otterlo, K.-H. Wagenblast, R. Fazio, and G. 
Schon, ibid. 48, 3316 (1993). 
[11] E. Granato and J.M. Kosterlitz, Phys. Rev. Lett. 65, 
1267 (1990). 

[12] C. Bruder, R. Fazio, A. Kampf, A. van Otterlo, and G. 
Schon, Phys. Scr. T 42, 159 (1992). 

[13] M.Y. Choi and S. Doniach, Phys. Rev. B 31, 4516 (1985). 

[14] The critical behavior of the system should not be affected 
by the choice of St; see, e.g., S. L. Sondhi, S. M. Girvin, 
J. P. Carini and D. Shahar, Rev. Mod. Phys. 69, 315 
(1997). Here all the dynamics of the system occurs over 
the time scale h/^8EjE (or H/yJWjFh for Co = 0), 
making the choice proper. 

[15] R. Savit, Rev. Mod. Phys. 52, 453 (1980); P.R. Thomas 
and M. Stone, Nucl. Phys. B 144, 513 (1978). 

[16] M.-S. Choi, J. Yi, M.Y. Choi, J. Choi, and S.-I. Lee, 
Phys. Rev. B 57, R716 (1998). 

[17] See, e.g., G. Blatter, M.V. Feigel'man, V.B. Geshken- 
bein, A.I. Larkin, and V.M. Vinokur, Rev. Mod. Phys. 
66, 1125 (1994), and references therein. 

[18] The upper critical field H C 2 for a type-II superconductor 
is irrelevant in our consideration here; ideally, the system 
corresponds to the limit H C 2 — > oo. 

[19] S. E. Korshunov, Europhys. Lett. 11, 757 (1990), and 
references therein. 

[20] P. Delsing, CD. Chen, D.B. Haviland, T. Bergsten, and 
T. Claeson, Lecture note in the Euroschool on Supercon- 
ductivity in Network and Mesoscopic Structures (unpub- 
lished); A.D Zaikin and S.V. Panyukov, Czech. J. Phys. 
46(S2), 629 (1996). 




FIG. 1. Superconducting array of size 9x6 with the mag- 
netic frustration / = 1/3. The 3x3 superlattices are indi- 
cated by thick solid lines. The position of the ith lattice site 
is represented by R + a, where R and a denote the position 
of the superlattice and the relative position of the site inside 
the superlattice, respectively. 
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FIG. 2. Three cases in which the second-order matrix Qij 
has nonzero values. Only the single-charge states \ j) with j 
at positions marked by the empty circle, the filled circles, and 
the 'x' symbols, can be connected by two successive charge 
hoppings, to the state \i) with i at the center (o). 
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FIG. 3. Phase boundaries between the Mott insulating 
phase (below each curve) and the superconducting phase 
(above), computed from the consideration of single-charge 
excitations. Boundaries for various ratios of the junction ca- 
pacitance Ci to the self-capacitance Co are shown: Ci/Co = 
0.0001(D), 0.1(o), 0.2(A), and 1.0(O). It is observed that the 
superconducting region expands as the junction capacitance 
is increased. 
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FIG. 4. The charge-dipole-type states \k,l) giving nonzero 
elements of the second-order matrix (i,j\M\k,l) in case that 
\i,j) has the charge configuration of a positive charge at site i 
and a negative charge at j as shown in (a). The filled and the 
empty circles denote the positive and the negative charges, 
respectively, (i) shows an example of the states not included 
in (a)-(h). 
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FIG. 5. Phase boundaries between the Mott insulating 
phase (below each curve) and the superconducting phase 
(above) of the charge-dipole type excitation. The upper curve 
is for Ci/C = 10000 while the lower one is for Ci/C = 10. 
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FIG. 7. Schematic diagram of the phase boundaries (a) for 
Co =fc and (b) for Co = 0. The dashed lines are valid only 
for continuum systems. 
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FIG. 6. Phase boundaries in the plane of Ej/Ei and / 
for Ci/Co — 100 between the Mott insulating phase (below 
the curve) and the superconducting phase (above) computed 
from the consideration of single-charge excitations. The ob- 
tained critical values (Ej/Ei) c w 0.657 for / = and 0.915 
for / = 1/2 are in good agreement with the experimental 
results. 
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